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Interferometric observations of the radio source pair 3C 84 and OE 400 have been 
made on the 21-km baseline between DSS 13 and DSS 15 to explore the angular 
navigation potential of intracomplex connected element interferometry (CEI). The 
differential phase-delay observable formed from pairs of 3- minute scans exhibited a 
precision of 1 psec, while the actual scatter of the phase-delay residuals for eleven 
scans over the 90-minute observing session was about 10 psec, consistent with the 
expected few-miilimeter fluctuations in the wet tropospheric path delay. Fitting 
for the position of OE 400 relative to 3C 84 yielded an error ellipse with a semi- 
minor axis of 60 nrad. Given the short data arc in this experiment, the orthogonal 
direction in the plane of the sky is not well determined; however , a second base- 
line or a data arc spanning a larger fraction of the source mutual visibility window 
could provide simultaneous determination of both right ascension and declination . 
Examination of the phase-delay residuals supports the accuracy of the cycle ambi- 
guity resolution. However , reliable phase ambiguity resolution will pose the most 
significant challenge to routine use of CEI for spacecraft tracking , particularly when 
the a priori spacecraft source position is not well known. Several approaches for 
ambiguity resolution are briefly outlined. 


I. Introduction 

Connected element interferometry (CEI) can provide 
accurate angular tracking from short intracomplex base- 
lines by making use of the very precise phase-delay data 
type. A number of operational advantages result from 
being able to form the angular tracking observable from 
data collected within a single Deep Space Communica- 
tions Complex (DSCC), including the potential for real- 
time data processing and fringe verification, as well as 
near-real-time delivery of tracking observables to the orbit 
determination navigation software. Some of the motiva- 


tions for CEI are discussed in [1], Previous intracomplex 
observations at Goldstone on the 6-km baseline between 
DSS 13 and DSS 12 have been reported [1,2]. 

CEI could provide an efficient angular tracking capa- 
bility during mission phases where the full precision of 
intercontinental very long baseline interferometry (VLBI) 
is not required. A 50-100 nrad intracomplex angular 
tracking capability on baselines of 20 km would represent 
an improvement over the current few-hundred-nrad angu- 
lar accuracy of Doppler tracking. In addition, the CEI 
observable would not require long, continuous data arcs 
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and would not suffer from the well-known singularity of 
the Doppler data type near zero declination. Increasing 
the CEI baseline lengths to 100 km could further improve 
angular accuracy to the 10-20 nrad level. Intercontinental 
VLBI would still provide the highest angular capabilities 
for specific mission requirements. Current VLBI accuracy 
is at the 30-50 nrad level for the operational Block I VLBI 
system, 1 while future systems may deliver nanoradian- 
level accuracy relative to the quasar frame [3] and 10-prad 
accuracy for relative angular tracking of two or more space- 
craft within the same primary antenna beam [4]. 

In this article, the first set of differential phase-delay 
observations on the 21-km baseline between DSS 13 and 
DSS 15 is presented. This is the longest baseline currently 
available within any of the DSCCs. DSS 13 is the Deep 
Space Network (DSN) Advanced Systems Program’s 26-m 
research and development antenna, while DSS 15 is a high- 
efficiency, 34-m antenna used primarily for operational 
DSN spacecraft tracking. Differential spacecraft-quasar 
observations have been simulated by observing a pair of 
quasars, 3C 84 and OE 400, separated by 6 deg in the 
plane of the sky. After describing the observing schedule 
and the data acquisition system, the observable formation 
will be reviewed. The resulting differential phase-delay 
residuals will be examined in light of previously reported 
calculations which suggest that tropospheric fluctuations 
should dominate the CEI differential phase observable [5]. 

Weighting the data to reflect the observed scatter, the 
differential phase-delay residuals will be used to solve for 
the relative positions of the two quasars. The resulting 
source position error ellipse characterizes the angular in- 
formation content of this intracomplex CEI pass. Since 
the a priori relative source positions are well known and 
based on many years of high-accuracy VLBI observations, 
the size of the estimated relative position shift also serves 
as a consistency check on the CEI results. 

Finally, potential problems relating to reliable ambigu- 
ity resolution will be discussed, along with a brief outline 
of several alternative approaches to ambiguity resolution. 

II. Experiment Description 

Previous intracomplex phase-delay observations at 
Goldstone have been made on the 6-km baseline between 


1 J. B. Thomas, “An Error Analysis for Galileo Angular Position 
Measurements with the Block I ADOR System,” JPL Engineering 
Memorandum 335-26 (internal document), Jet Propulsion Labora- 
tory, Pasadena, California, November 11, 1981. 


DSS 12 and DSS 13, using a fiber-optic link between those 
stations to operate them coherently [1,2]. These obser- 
vations were encouraging in that phase connection was 
reliably achieved and the final phase-delay residuals were 
just a few millimeters, but the extremely short baseline 
limited the angular resolution of these observations to 
many hundreds of nrad. When the fiber-optic link at 
Goldstone was extended to SPC-10, where both DSS 14 
and DSS 15 are located, a 21-km baseline between DSS 13 
and the SPC-10 antennas became available. This longer 
baseline offers the potential for much greater angular accu- 
racy. With the goals of understanding the limiting phase- 
delay errors, demonstrating phase ambiguity resolution, 
and quantifying the angular accuracy on this 21-km base- 
line, an experiment was planned and scheduled. 

Interferometric data were acquired on May 27, 1988, 
from 2130-2300 UT, at DSS 13 and DSS 15 in the Gold- 
stone DSCC. Observations were made at both 2.3 GHz (S- 
band) and 8.4 GHz (X-band) to allow calibration of the 
effects of charged particles. A hydrogen maser at SPC-10 
provided the frequency reference for DSS 15; the analog 
fiber-optic link was used to transfer this frequency refer- 
ence to DSS 13, allowing the two stations to be operated 
coherently. Previously reported tests of the fiber-optic link 
indicate that the link stability A vjv is well below 10~ 14 
for averaging times of several hundred seconds [6]. 

The Block 0 VLBI data acquisition system was used to 
single-bit quantize, sample, format, and record the incom- 
ing signals at each station [7]. The Block 0 system records 
a single 2-MHz baseband signal, sampled with 1-bit quan- 
tization at the Nyquist rate, resulting in a 4-Mbit/sec data 
rate. This 2-MHz bandwidth was time-multiplexed over 
six separate sky frequencies, three at S-band and three at 
X-band. Each S-band channel was observed for 0.4 sec, 
while each X-band channel was observed for 1.6 sec, dur- 
ing each 6-second multiplexing cycle. (The longer X-band 
dwells were chosen because the X-band observable carries 
a greater weight in the final S/X linear combination.) 

A. Frequency Constraints on Short Baselines 

Table 1 shows the sky frequencies and channel dwell 
times used in this experiment. Note the 200-IIz offset be- 
tween the local oscillator (LO) frequencies at the two sta- 
tions. This is required due to the small differential Doppler 
shift between stations, and can be understood as follows: 
The sidereal rotation of the Earth induces a Doppler shift, 
for observations in the direction of a given radio source, at 
each station of an interferometer. For short baseline ob- 
servations, the differential Doppler shift v D between the 
two stations becomes very small, typically less than 10 Hz 
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at 2.3 GHz on a 21-km baseline. Adding a local oscilla- 
tor offset vqff ensures that when the signals from the 
two stations are cross-correlated (i.e., multiplied), the re- 
sulting fringe frequency, up = ^ + i^oFF) will be well 
separated from dc. When the fringes are subsequently 
counter-rotated, or “stopped,” by multiplying them with 
a model sinusoid with frequency vp « up, both sum and 
difference frequencies are generated, at roughly 2 up and 0, 
respectively. In the subsequent coherent integration, the 
high-frequency term will be attenuated if the integration 
time r 1/(2^f), leaving only the stopped fringe phase, 
as desired. With the coherent integration time of 0.2 sec 
used at the correlator, the 200-Hz LO shift ensures that 
this is the case. 

Phase calibration tones were injected at each station 
at both S-band and X-band frequencies to calibrate tem- 
poral fluctuations in the instrumental phase. Here again, 
some modifications to standard VLBI procedure were re- 
quired due to the short baselines. If the same calibra- 
tion tone frequencies are used at both stations, the cross- 
correlation of these phase calibration tones will produce a 
signal with very nearly the same frequency as the actual 
quasar fringes, even when an LO offset is used, due to the 
very small differential Doppler shift between stations on a 
short baseline. For this reason, different calibration tone 
frequencies were used at each station in this experiment. 
The DSN phase calibration system produces calibration 
tones at integral multiples of b/N MHz, where N is user- 
selectable in the range 5-99. A value of N = 7 was used 
for DSS 13, and TV = 8 for DSS 15. This choice placed at 
least three calibration tones in each 2-MHz channel. All 
tones in each channel were phase-tracked and used dur- 
ing the correlation process to calibrate instrumental phase 
errors. 

B. Observations and Initial Processing 

The experiment observing schedule consisted of re- 
peated differential observations of the two radio sources 
3C 84 and OE 400. These sources are separated by about 
6 degrees in the plane of the sky. Table 2 summarizes 
the a priori J2000 positions of the two sources, based on 
many years of DSN VLBI observations. Each source was 
observed for 180 sec, with a slew time of 30-45 sec be- 
tween observations. Both sources were setting during the 
period of observation; 3C 84 varied in elevation from 57° 
down to 42°, while OE 400 varied from 53° down to 40°. 
(On the DSS 13-DSS 15 baseline, the elevation angles at 
the two stations will always agree to 0.2 deg or better.) 
Table 3 lists the sequence of observations. One observation 
pair was lost due to a brief failure of the DSS 13 antenna 
pointing computer at 215945 GMT. All other scans were 
successful. 


The recorded data, on Block 0 VLBI videocassettes, 
were processed at the JPL/Caltech Block II VLBI correla- 
tor [8]. The output of this initial processing is sine and co- 
sine correlation sums for each of the six channels at a rate 
of one point per channel every six seconds, corresponding 
to the time-multiplexing cycle period. These correlation 
sums were then processed with the REVERT/PII ASOR 
software package to fit for fringe amplitude and phase for 
each channel. Fringes were obtained for all of the obser- 
vations for which data were obtained. 


With three frequency channels recorded at S-band and 
three at X-band, it is possible to form group delay ob- 
servables using the bandwidth synthesis (BWS) technique 
[9]. Although not a goal of this experiment, the BWS ob- 
servable was formed as a check on the quality of the data. 
The BWS observable is formed by calculating the slope of 
phase versus frequency over a given spanned bandwidth: 


tbws — 
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For the channel frequencies used in this experiment, the 
maximum spanned bandwidth v\ —V 2 was 40 MHz at both 
S-band and X-band. The formal error on the X-band BWS 
delay observable was about 100 psec for 3C 84, and over 
200 psec for OE 400, with somewhat higher errors at S- 
band due to the shorter dwell times. These errors corre- 
spond to path-delay errors on the order of 5 cm, or angular 
errors of several /irad on a 21-km baseline. 


III. The Phase Observable and Ambiguity 
Resolution 

To achieve angular accuracies of better than 100 nrad, 
the much more precise phase data type is required. The 
geometric component of the interferometric phase can best 
be thought of as a measure of the geometric delay in units 
of the observing wavelength. To make use of the interfer- 
ometric phase, however, one must be able to resolve the 
integer cycle ambiguity associated with that data type. 
Further complicating the ambiguity resolution is the con- 
tribution of nongeometric errors to the phase, such as the 
unknown LO phase offset between stations, propagation 
media delays, and unknown instrumental delays and dis- 
persions. 

In principle, the BWS delay residuals contain infor- 
mation about errors in the correlator delay model that 
could be used to resolve the phase observable on a scan- 
by-scan basis. However, the formal errors on the BWS 
delay residuals in this experiment were on the order of an 
RF cycle (120 psec at X-band) or larger, and thus could 
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not be used to aid in ambiguity resolution. (It should be 
pointed out, however, that a wider BWS spanned band- 
width coupled with a larger total recorded bandwidth 
could provide sufficient BWS precision to enable cycle am- 
biguity resolution for each scan based on the BWS resid- 
ual. This more robust approach to ambiguity resolution 
will be attempted in an experiment planned for early 1990 
on this same baseline, using the wider bandwidth Block II 
VLBI system, which is capable of providing 10-psec BWS 
precision.) 

Instead, for this experiment it was necessary to deter- 
mine an a priori delay model which was accurate enough to 
resolve the cycle ambiguity directly. Using the MASTER- 
FIT VLBI parameter estimation software [10] and the best 
a priori models for station locations and source positions, 
an (ambiguous) residual phase observable was calculated 
for each quasar observation. Consider observations of two 
sources, A and D ) and let <j>A represent the total phase 
observable for source A, Na the unknown integer cycle 
ambiguity, wrf the RF observing frequency, and ta the 
a priori model delay. Then the phase residual 4>a can be 
written 

<t>A ~ <f>A + 2ttNa — UrfTa 

A similar phase residual can be formed for the subsequent 
observation of source B: 

4>b = <f>B ± 27 tNb — wrftb 

If the a priori model were sufficiently accurate, the integer 
cycle ambiguities N a and Nr could each be determined 
by the requirement that the phase residuals be near zero. 
However, the above-mentioned nongeometric model uncer- 
tainties are typically on the order of an RF cycle or more, 
generally preventing this. 

By differencing the two phase observations, many of 
the model uncertainties are canceled or reduced, thereby 
allowing the determination of the relative cycle ambiguity 
Na — Nr. Forming a differential observable yields 

= (<1>a ~ $b) + 2tt(Na - N b ) - urf(ta - r B ) 

= A<(>a-b + ZttANa-b — wrfAta-b 

Many potential error sources in the delay models f are 
reduced or eliminated by this differencing. The unknown 
offset of the LO phases between the tw r o stations, which 
manifests itself as a clock offset for the single-source ob- 
servable, cancels completely in the differential observable. 
Geometric errors are also largely reduced by differencing 
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the phase delays for angularly close sources. A baseline 
uncertainty of size A B can cause an error of up to AB/c 
in the delay model for an individual source, where c is the 
speed of light. However, when two sources with a small 
angular separation of A 8 are differenced, much of this ge- 
ometric error cancels. The resulting error in the differ- 
enced model delay will be bounded by ABAQfc, where 
A# is expressed in radians. For example, a 6-deg angular 
source separation will cause geometric model errors to be 
reduced by about a factor of 10 in a differential observable. 
Thus even a 2-cm baseline error, representing a full half- 
cycle error for X-band observations, would be reduced to a 
2-mm maximum path-delay error for a differential observ- 
able between sources separated by 6 deg, and would not 
complicate the relative phase connection. 

Similarly, propagation media errors are greatly reduced 
by differencing observables for angularly close sources. 
Any uncertainty in the overall tropospheric delay at each 
station is greatly attenuated due to several factors: The 
short baseline causes the total zenith delays at the two sta- 
tions to be highly correlated; the short baseline also causes 
the two stations to observe a source at nearly the same ele- 
vation angle; and the angular proximity of the sources also 
causes the elevation angles for the two sources to be nearly 
the same at each station. The remaining error is predomi- 
nantly due to small-scale temporal and spatial fluctuations 
in the tropospheric delay, on the time scale of the time 
between scans, and on the spatial scale of the angular dis- 
tance between sources, projected to the tropospheric scale 
height of several kilometers. 

The relative cycle ambiguity AN is given by 

AAT = nint p FA ^- A * ) 

where nint is the “nearest integer” function, and where 
the subscript A — B has been dropped. The reliability of 
the phase ambiguity determination is reflected in the size 
of the final phase residuals. By the above choice of AN, 
the residuals will lie between ±1/2 cycle; for the phase 
connection to be deemed reliable, however, the residuals 
should be distributed in a peak around zero and be well 
separated from ±1/2 cycle. 

For this condition to be satisfied, two criteria must be 
met: both the precision of the differential phase observ- 
able A<j) and the error in the a priori differential model 
delay At must be well below 1/2 cycle of phase. The first 
criterion is usually satisfied just by the requirement that 
the sources be detected with good SNR since the statis- 
tical phase error, expressed in radians, is roughly 1/SNR. 


i 



The second criterion is more difficult to ensure. Given 
known uncertainties in station location and source posi- 
tion, the geometric component of A f can be calculated re- 
liably. However, contributions of unmodeled errors, such 
as stochastic troposphere fluctuations or antenna defor- 
mation, are more difficult, to estimate. The approach in 
this article is to simply determine the phase connection 
using the best a priori information, and then evaluate the 
reliability of the phase connection by examining the dis- 
tribution of phase residuals A <f>. 

Figure 1 shows histograms of the S-band and X-band 
phase residuals after ambiguity resolution. The residuals 
show a very clear clustering about zero, with a root-mean- 
square (rms) value below l/10th of a cycle for both bands. 
Quantitatively determining the reliability of the ambiguity 
resolution requires some assumption about the underlying 
probability distribution of the phase residual error sources. 
If it were assumed that the phase residual error was due 
to a Gaussian error source with a standard deviation of 
0.1 cycle, then it would be possible to calculate the like- 
lihood that all of the integer cycle ambiguities had been 
correctly resolved. The probability that a single observa- 
tion was correctly resolved would simply be the probability 
that the observation’s phase error was less than 0.5 cycle. 
For the case of (T = 0.1 cycle, this is just the integrated 
Gaussian probability distribution out to ±5 <7, equal to 
0.9999994. For the eleven observations, this gives a cu- 
mulative probability of 0.999993 that all eleven observa- 
tions have been correctly resolved, if it is assumed that 
the observations are uncorrelated. Unfortunately, the er- 
ror sources contributing to the observed phase residuals 
most likely have significant nonGaussian tails that would 
drastically change this conclusion, increasing the probabil- 
ity of an incorrect integer cycle determination. Obtaining 
a better characterization of the probability distribution of 
phase errors, and thus determining more realistically the 
reliability of cycle ambiguity resolution, is a primary mo- 
tivation for collecting a much larger database of observa- 
tions on this intracomplex baseline. At this time, all that 
can be said is that the distribution of phase residuals af- 
ter ambiguity resolution strongly suggests that the phase 
ambiguity resolution was successful for this experiment. 

In the final stage of observable formation, the S-band 
and X-band phase residuals are linearly combined to form 
an S/X residual phase-delay observable, free from the dis- 
persive effects of charged particles: 
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Figure 2 shows the S/X-corrected differential phase- 
delay residuals as a function of time for the eleven source 
pair observations. The rms value of the residuals is 10 psec, 
or about 3 mm of path delay. The formal statistical error 
on each point is only about 1 psec; the observed scatter 
is expected to be due primarily to stochastic spatial and 
temporal fluctuations of the wet tropospheric path delay 
at each station. To account for these tropospheric fluc- 
tuations, the actual error bars shown on the data points 
have been inflated to yield a reduced \ 2 equal to 1, as 
described in Section IV. A previous article [5] calculated 
the size of such fluctuations and their impact on differen- 
tial phase-delay observations in CEI, based on numerical 
integrations of a Kolmogorov turbulence model for atmo- 
spheric fluctuations [11]. Using the model described in 
that article and the relevant parameters of these observa- 
tions (namely, a 21-km baseline, 6-deg source separation, 
3-minute scan and 40-sec slew time, and the range of el- 
evation angles occurring during this experiment) the ex- 
pected scatter in the differential phase delay is calculated 
to range from 3.5 to 4.2 mm over the observation period, 
increasing as the sources set in elevation. The model has 
some uncertainty due to daily variations in the wind speed 
and scale height of the wet troposphere. Nonetheless, this 
excellent level of agreement supports the hypothesis that 
tropospheric effects are indeed a dominant error source for 
the CEI differential phase observable. 

IV. Relative Source Position Accuracy 

The most direct way of demonstrating the angular ac- 
curacy of these observations is to solve for the relative an- 
gular positions of the two sources, just as one would solve 
for a spacecraft position relative to a reference quasar. Do- 
ing this, one obtains an estimated position correction for 
one of the quasars, along with a covariance describing the 
uncertainty in that estimated position. The a priori source 
positions of both sources in this experiment are known to 
an accuracy of about 10 nrad, based on many years of 
observations on long baselines. Thus a consistency check 
on the estimated position is that it agree, to within the 
estimated covariance, with the a priori position. 

(There may be concern that structure in the radio 
source brightness distribution could cause the apparent 
source position to differ on long versus short baselines. 
Although source structure issues will not be addressed 
in this article, future CEI observations of a number of 
source pairs would provide a data set for quantifying 
this effect. If long- and short-baseline apparent source 
positions differed significantly, it would be necessary to 
develop a distinct CEI source position catalog based solely 
on short-baseline observations.) 
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The MASTERFIT VLBI software was again used, in 
this case to estimate a position correction to the source 
OE 400. As mentioned earlier, the formal statistical er- 
rors on the phase-delay data points are only about 1 psec, 
while the actual scatter of 10 psec is consistent with the 
expected level of tropospheric fluctuations. To take this 
into account, an error contribution of about 10 psec was 
added in quadrature to each differential observation’s for- 
mal error to obtain the final data weight used in the source 
position parameter estimation. The size of this supplemen- 
tal error contribution was determined by requiring the re- 
duced x 2 °f the final fit to equal one. In this way, the 
actual level of fluctuations observed in the data can be in- 
corporated into the data weights and thereby into the un- 
certainty in the final source position determination. Given 
that atmospheric dynamics can vary significantly from day 
to day, this approach is probably more robust and reliable 
than trying to specify the data weights based solely on 
some model of tropospheric fluctuations. Nonetheless, as 
described earlier, this supplemental phase-delay error con- 
tribution, determined empirically, agrees fairly well with 
the predictions of a statistical model of tropospheric fluc- 
tuations. 

Using these empirically determined data weights, the 
eleven differential phase-delay residuals served as input to 
a weighted least-squares adjustment of the right ascension 
and declination of OE 400. No other parameters were 
estimated in this process. The OE 400 source position 
was essentially unconstrained: the a priori right ascension 
and declination uncertainties were set at 1 radian. Fig- 
ure 3 shows the resulting \-<r error ellipse for the estimated 
OE 400 position correction. The vertical axis is the dec- 
lination shift A <5, while the horizontal axis represents the 
right ascension shift expressed as an arclength: Ao x cos 6. 
The error ellipse is very elongated, with a semi-major axis 
of about 1165 nrad, due to the short data arc. However, in 
the direction of the average baseline projection during the 
experiment, the semi-minor axis is 60 nrad. In addition, 
the l-<7 error ellipse is consistent with the a priori source 
position, which in this experiment was well known. 

These results represent the first time that angular mea- 
surements from a single DSN complex have provided ac- 
curacies below 100 nrad. The high eccentricity of e = 0.90 
for the error ellipse derived from this data set reflects the 
short span of data. Each single phase-delay observation 
only contains information for one direction in the plane 
of the sky, namely, along the projection of the interfer- 
ometer baseline in the radio source direction. Over the 
90-minute duration of this experiment, the baseline pro- 
jection in the direction toward OE 400 rotates only a small 
amount, about 10°. Because the experiment duration was 


so short, the full range of baseline orientation was not sam- 
pled, and thus the resulting source position will have a 
much better determination in one direction — namely, the 
projected baseline direction at the central epoch of the ex- 
periment, with a much larger uncertainty in the orthogonal 
direction. 

To illustrate this, Fig. 4 shows the east-west and north- 
south components of the projected DSS 13-DSS 15 base- 
line in the direction of OE 400, over the full 15. 5-hour mu- 
tual visibility window of this northern declination source. 
(This can be thought of as the baseline length and orien- 
tation, as viewed from the quasar.) The 90-minute ob- 
servation period represented by the current data set is 
indicated by the heavy line. The limited range of base- 
line orientation during the 90-minute observation period 
underlies the large eccentricity of the estimated OE 400 
position error ellipse in this experiment. But over the en- 
tire mutual visibility window, the baseline rotates through 
well over 90°, and hence CEI observations spanning this 
entire period would provide strong solutions for both com- 
ponents of sky position. Previous analyses 2,3 have shown 
that CEI observations at Goldstone on the DSS 13-DSS 15 
baseline, collected throughout the mutual visibility period, 
can provide good determination of both a and 6 for north- 
ern declination sources. 

Adding a second baseline, orthogonal to the DSS 13- 
DSS 15 baseline, would also enable good determination of 
both components of sky position. This solution has the 
added advantage that both components could be deter- 
mined simultaneously from observations taken at a single 
epoch, instead of combining observations that are sepa- 
rated by many hours—an important advantage in the case 
of tracking a spacecraft with a poorly determined orbit. 

V. Reliability of Ambiguity Resolution 

Successful cycle ambiguity resolution is the key to un- 
locking the high precision of the phase data type. Much 
more work needs to be done to understand the reliability 
of ambiguity resolution as a function of baseline length, 
angular source separation, elevation angle, and temporal 
scan separation. Nonetheless, one encouraging conclusion 


2 M. II. Finger and C. D. Edwards, “Relative CEI Navigation Per- 
formance of Goldstone Intracomplex Baselines," JPL Interoffice 
Memorandum 335.3-88-116 (internal document), Jet Propulsion 
Laboratory, Pasadena, California, October 20, 1988. 

3 S. W. Thurman, “Information Content of a Single Pass of Phase- 
Delay Data From a Short Baseline Connected Element Interfer- 
ometer," JPL Engineering Memorandum 314-479 (internal docu- 
ment), Jet Propulsion Laboratory, Pasadena, California, December 
13, 1989. 
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of this analysis is that differential phase-delay errors do 
not seem to grow linearly with baseline length. Limited 
data on a 6-km, a 21-km, and a 253-km baseline are now 
available, and in each case, the final delay residual scat- 
ter has been at about the 10-psec level. This is consis- 
tent with the hypothesis that the scatter is dominated at 
each site by rapid temporal fluctuations, which are un- 
correlated between stations on the short time scale of a 
differential observation, and therefore are independent of 
baseline length. In each of these experiments, the a pri- 
ori source positions were well known; thus the ambiguity 
resolution was only limited by the size of unmodeled com- 
ponents of the differential delay, particularly these rapid 
tropospheric fluctuations. At S-band and X-band, this 10- 
psec level of unmodeled phase-delay uncertainty does not 
prevent ambiguity resolution. However, at 32 GHz (Ka- 
band), the delay ambiguity is only about 30 psec. These 
same errors would then represent a much larger fraction 
of an RF cycle, seriously complicating phase connection. 

Ambiguity resolution also becomes more difficult when 
the a priori source position of one of the sources is less 
well known. In the “direct” ambiguity resolution strategy 
used in this article, in which the integer cycle ambiguity is 
determined solely on the a priori delay model, one clearly 
must know the source position in the projected baseline 
direction to better than half of a fringe spacing. On a 
21-km baseline, the minimum fringe spacing is 6.2 /irad 
at 2.3 GHz and 1.7 //rad at 8.4 GHz. Doppler tracking 
can usually determine a spacecraft ephemeris with suffi- 
cient accuracy to satisfy this constraint during periods of 
cruise. However, at encounter or during other periods of 
high spacecraft dynamics, this constraint may be more dif- 
ficult to satisfy. 

Several approaches can be used to improve the relia- 
bility of ambiguity resolution in the case that the a pri- 
ori model delay error is too large to allow “direct” am- 
biguity resolution. All of them involve starting with a 
lower accuracy data type, but one whose cycle ambiguity 
can be reliably resolved. The information in this lower 
accuracy observable is then used to assist in resolving 
the cycle ambiguity associated with the higher accuracy 
data type. Some of these approaches have been discussed 
previously and are outlined below. 4 What follows is a 
sketch of possible techniques; more work is required to fully 
analyze and explore how each of them would perform in 
an operational CEI system. 


4 S. Thurman, “Continuation of the Connected Element Interfer- 
ometry System Study,” JPL Interoffice Memorandum 314.5-1367 
(internal document), Jet Propulsion Laboratory, Pasadena, Cali- 
fornia, October 20, 1988. 


A. BWS Aiding 

As mentioned previously, with sufficient spanned band- 
width, BWS group delay residuals can be determined with 
formal errors of less than a cycle of RF phase. In that 
case, if dispersive errors can be accurately calibrated and 
removed, the group delay residual can be used to correct 
the delay model and aid in ambiguity resolution. How 
precise must the BWS observable be to aid in ambiguity 
resolution? Consider a system with two BWS channels, 
separated in frequency by a spanned bandwidth of vbws , 
and with a formal phase error of cycles. Requiring that 
the BWS observable precision be less than l/6th of an RF 
cycle, to ensure reliable ambiguity resolution, yields: 

\/2<7 (p 1 

vbws Svrf 

With a 400- MIIz spanned bandwidth and a 5-mcyc phase 
error in each channel, it would then be possible to resolve 
RF cycles at X-band. It is important to keep in mind that 
this argument has ignored the effects of dispersive errors 
due to the ionosphere, solar plasma, or instrumentation, 
which contribute differently to the phase and group delay 
data types. 

For use in spacecraft tracking, this technique of BWS- 
aided ambiguity resolution will require large spanned 
band widths on the spacecraft. However, the current X- 
band deep space downlink frequency allocation is only 
40 MHz wide, probably insufficient for resolving the X- 
band or even the S-band phase ambiguity. At Ka-band, 
however, a full 500-MHz downlink allocation exists and 
could aid ambiguity resolution. For example, with a space- 
craft incorporating a dual-band X/Ka downlink, and with 
widely spaced VLBI tones at Ka-band, the accurate Ka- 
band BWS observable could be used to resolve the X-band 
phase observable. Such a scenario has been proposed for 
same-beamwidth interferometric tracking of two spacecraft 
at Mars with simultaneous X/Ka downlinks [12]. 

B. Band-to-Band Aiding 

Just as one can “step” from the BWS observable, with 
its spanned bandwidth vbws , u P to the phase observable 
at frequency i / rf , one can also imagine using the phase ob- 
servable at a lower frequency band to aid the phase ambi- 
guity resolution at a higher frequency band, e.g., using an 
S-band phase residual to resolve the X-band phase ambigu- 
ity. Again, the key to making this work is the requirement 
that dispersive errors represent much less than a cycle of 
phase at the higher frequency. Instrumental phase cali- 
brations can probably reduce instrumental dispersion to a 
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few degrees of phase, in which case they would not pose 
an obstacle to stepping from one band to the next. In that 
case, the biggest dispersive error will be due to the effects 
of charged particles in the ionosphere and solar plasma. 
Since charged-particle delays scale as 1/r/ 2 , whereas cycle 
ambiguities scale as 1/T, charged particles will represent a 
four times more serious error in stepping from S-band to 
X-band than in stepping from X-band to Ka-band. 

The current understanding of ionospheric errors and, in 
particular, spatial and temporal fluctuations in ionospheric 
delays is not as well developed as the corresponding un- 
derstanding of tropospheric errors. Empirically, however, 
the concept of S-to-X-band aiding in the current data set 
can be tested by asking the question: if the final S-band 
phase delays had been used instead of the a priori delay 
model to resolve the X-band phase, would the same inte- 
ger cycle ambiguities have been determined? Examining 
the S-band phase delay residuals, it was found that ten of 
the eleven X-band observables would have been resolved 
identically. However, one observation (the outlier in the 
S-band histogram of Fig. 1) would have had a 1-cycle shift 
in the final X-band observable. Presumably the “direct” 
solution in this case is the correct one, not the “S-to-X- 
aided” solution, based on the final distribution of X-band 
residuals for the repeated observations of the source pair. 
It is not known whether the S-band outlier was caused 
by an uncalibrated instrumental fluctuation or by actual 
fluctuations in the ionosphere above Goldstone. 

C. Baseline Aiding 

Going to shorter baselines increases the fringe spacing 
on the sky, and thus relaxes the a priori source position 
requirements for phase connection (at the cost of lower 
angular precision). One can imagine a series of progres- 
sively longer baselines, where the phase residual on each 
baseline would be used to aid the ambiguity resolution of 
the next longer baseline. If the a priori position knowl- 
edge is sufficient to resolve the phase ambiguity on the 
shortest baseline, and if each baseline’s phase-delay resid- 
ual is sufficiently accurate to resolve the ambiguity on the 
next longer baseline, then ultimately the phase ambiguity 
on the longest baseline can be resolved. This technique 
works best if the baselines are exactly parallel, since each 
baseline only measures instantaneously one component of 
the source position. 

Consider a linear array of five antennas with coordi- 
nates of 0, 1, 5, 25, and 125 km along some line. At 
X-band, the 1-km baseline formed by the first two anten- 
nas would correspond to a very large minimum fringe spac- 
ing of 40 //rad. Phase connection on this baseline would 


require only very crude a priori source position informa- 
tion. If the phase observable on this 1-km baseline were 
sufficient to determine the source position (along the pro- 
jected baseline direction) to a l-cr uncertainty of 1 / 30th 
of a fringe spacing, that improved source position would 
enable reliable ambiguity resolution on the 5-km baseline, 
i.e., the 3-<r source position error would contribute less 
than half a cycle of phase on the 5-km baseline. This pro- 
cess would continue until the 125- km baseline was reached. 
A 1 / 3 0 1 h of a cycle error on the 125-km baseline would 
correspond to a 10-nrad angular position determination at 
X-band. 

Each of these techniques has inherent strengths and 
weaknesses in terms of requirements of the CEI system. 
The BWS or band-to-band techniques could help to resolve 
ambiguities at Ka-band in the likely event that stochastic 
tropospheric fluctuations prevent direct Ka-band ambigu- 
ity resolution. All three techniques could help in the case 
of tracking an object with very poor a priori angular posi- 
tion information. Perhaps the most robust solution would 
utilize several of these approaches at once: for instance, 
an array of baselines receiving simultaneous X/Ka down- 
links, with widely spaced Ka-band tones for generating a 
high precision group delay data type. 

VI. Summary 

A 90-minute observation pass on the DSS 13-DSS 15 
baseline has yielded eleven differential phase-delay observ- 
ables with formal errors of I psec and an rrns scatter of 
10 psec, roughly the level predicted by statistical mod- 
els of tropospheric fluctuations. These data determined 
the relative source positions to an accuracy of 60 nrad in 
one component of sky position. This is the highest an- 
gular accuracy achieved to date from observations within 
a single DSN tracking complex. The other component of 
source position was not well determined, due to the short 
observing period. However, similar observations spanning 
the mutual visibility window for this source would provide 
comparable accuracy for both components of sky position. 
Placing a quantitative value on the probability of a cy- 
cle error is difficult without a better understanding of the 
probability distributions of the underlying error sources, 
although examination of the distribution of phase residu- 
als suggests that the cycle ambiguity resolution was suc- 
cessful. 

Developing techniques to ensure highly reliable ambi- 
guity resolution is a prerequisite to incorporating CEI 
into operational spacecraft tracking. At Ka-band, the 
few-millimeter level of stochastic tropospheric fluctuations 
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observed in this experiment would seriously complicate 
ambiguity resolution. In addition, large uncertainties in 
a priori source position would also hamper ambiguity res- 
olution. Several techniques have been outlined for using 


less accurate, but more easily resolved, data types to aid in 
resolving a higher accuracy observable. Data collected in 
the last few months and over the coming year should allow 
some of these techniques to be tested and evaluated. 
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Table 1. Frequency configuration 


Channel 

DSS 13 frequency, 
MHz 

DSS 15 frequency, 
MHz 

Sideband, 

upper/lower 

Dwell time, 
sec 

1 

2304.9999 

2305.0001 

LSB 

0.4 

2 

2264.9999 

2265.0001 

USB 

0.4 

3 

2297.9999 

2298.0001 

LSB 

0.4 

4 

8439.9999 

8440.0001 

LSB 

1.6 

5 

8399.9999 

8400.0001 

USB 

1.6 

6 

8432.9999 

8433.0001 

LSB 

1.6 


Table 2. The a priori J2000 source positions 

Source 

Right ascension, 
h m s 

Declination, 
o / // 

3C 84 

03 19 48.160328 

41 30 42.10581 

OE 400 

03 03 35.242181 

47 16 16.27738 


Table 3. Observation sequence 


Observation no. 

Source name 

Start — stop, 
GMT 

Average elevation, 
deg 

1 

3C 84 

213100 - 213400 

56.7 

2 

OE 400 

213442 - 213742 

53.2 

3 

3C 84 

213824 - 214124 

55.3 

4 

OE 400 

214206 - 214506 

51.9 

5 

3C 84 

214542 - 214842 

54.9 

6 

OE 400 

214924 - 215224 

50.7 

7 

3C 84 

215300 - 215600 

52.6 

8 

OE 400 

215642 - 215942 

49.5 

9 a 

3C 84 

220018 - 220318 

51.2 

10 a 

OE 400 

220354 - 220654 

48.2 

11 

3C 84 

220730 - 221030 

49.9 

12 

OE 400 

221106 - 221406 

47.0 

13 

3C 84 

221442 - 221742 

48.5 

14 

OE 400 

221818 - 222118 

45.8 

15 

3C 84 

222154 - 222454 

47.2 

16 

OE 400 

222530 - 222830 

44.6 

17 

3C 84 

222906 - 223206 

45.9 

18 

OE 400 

223242 - 223542 

43.4 

19 

3C 84 

223618 - 223918 

44.5 

20 

OE 400 

223954 - 224254 

42.2 

21 

3C 84 

224330 - 224630 

43.2 

22 

OE 400 

224706 - 225006 

40.9 

23 

3C 84 

225042 - 225342 

41.9 

24 

OE 400 

225418 - 225718 

39.7 


These two observations were lost due to a malfunction of the antenna pointing computer at 
DSS 13. (Observations are grouped in pairs used to form the final differential observables.) 
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Fig. 1. S-band and X-band phase residuals after ambiguity reso- 
lution. The rms width of the distribution is 0.087 cycle at S-band 
and 0.093 cycle at X-band. 



Fig. 2. Final S/X-combined differential phase-delay residuals as 
a function of time over the 90-minute observation. 



Fig. 3. Estimated source position correction and error ellipse 
for OE 400, based on a weighted least-squares fit to the phase- 
delay residuals shown in Fig. 2. The orientation of the ellipse 
corresponds to the average projected baseline direction over the 
90-minute observation period. 



Fig. 4. The projection of the DSS 13-DSS 15 baseline on the plane 
of the sky, in the direction of the radio source OE 400, over the full 
mutual visibility window of over 15 hours. The portion of the vis- 
ibility window represented by the 90-minute observation reported 
here is shown with a heavier line; during this short period, the 
projected baseline only rotates through about 10 degrees, and 
so the resulting position determination is much stronger in one 
direction on the plane of the sky than in the other. 
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